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EXISTENCE OF GLOBAL WEAK SOLUTIONS FOR 3D 
DEGENERATE COMPRESSIBLE NAVIER-STOKES EQUATIONS 

ALEXIS F. VASSEUR AND CHENG YU 


Abstract. In this paper, we prove the existence of global weak solutions for 3D com¬ 
pressible Navier-Stokes equations with degenerate viscosity. The method is based on the 
Bresch and Desjardins entropy conservation [2] - The main contribution of this paper is 
to derive the Mellet-Vasseur type inequality [32] for the weak solutions, even if it is not 
verified by the first level of approximation. This provides existence of global solutions 
in time, for the compressible Navier-Stokes equations, for any 7 > 1 in two dimensional 
space and for 1 < 7 < 3 in three dimensional space, with large initial data possibly 
vanishing on the vacuum. This solves an open problem proposed by Lions in m 


1. Introduction 

The existence of global weak solutions of compressible Navier-Stokes equations with 
degenerate viscosity has been a long standing open problem. The objective of this current 
paper is to establish the existence of global weak solutions to the following 3D compressible 
Navier-Stokes equations: 

p t + div(pu) =0 

(pu)t + div(pu ®u)+VP- div(pOu) = 0, 

with initial data 

p\t=o = Po(x), pu | t=0 = m 0 (x). (1.2) 

where P = p 7 , 7 > 1, denotes the pressure, p is the density of fluid, u stands for the 
velocity of fluid, Bu = |[Vu +V 7 u] is the strain tensor. For the sake of simplicity we will 
consider the case of bounded domain with periodic boundary conditions, namely D = T 3 . 

In the case 7 = 2 in two dimensional space, this corresponds to the shallow water 
equations, where p(t, x ) stands for the height of the water at position x, and time t, and 
u(f, x ) is the 2D velocity at the same position, and same time. In this case, the physical 
viscosity was formally derived as in m (see Gent [H]). In this context, the global 
existence of weak solutions to equations (11.11) is proposed as an open problem by Lions in 
m ■ A careful derivation of the shallow water equations with the following viscosity term 

2 div(pBu) + 2 V(pdivu) 
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can be found in the recent work by Marche [28]. Bresch-Noble mm provided the mathe¬ 
matical derivation of viscous shallow-water equations with the above viscosity. However, 
this viscosity cannot be covered by the BD entropy. 

Compared with the incompressible flows, dealing with the vacuum is a very challenging 
problem in the study of the compressible flows. Kazhikhov and Shelukhin [25] established 
the first existence result on the compressible Navier-Stokes equations in one dimensional 
space. Due to the difficulty from the vacuum, the initial density should be bounded away 
from zero in their work. It has been extended by Serre [ 33 ] and Hoff [20] for the discon¬ 
tinuous initial data, and by Mellet-Vasseur [33] in the case of density dependent viscosity 
coefficient, see also in spherically symmetric case dm mi ns]. For the multidimensional 
case, Matsumura and Nishida mmm first established the global existence with the 
small initial data, and later by Hoff [211 122] [23] for discontinuous initial data. To re¬ 
move the difficulty from the vacuum, Lions in [27] introduced the concept of renormalized 
solutions to establish the global existence of weak solutions for 7 > | concerning large 
initial data that may vanish, and then Feireisl-Novotny-Petzeltova m and Feireisl da 
extended the existence results to 7 > |, and even to Navier-Stokes-Fourier system. In all 
above works, the viscosity coefficients were assumed to be fixed positive numbers. This is 
important to control the gradient of the velocity, in the context of solutions close to an 
equilibrium, a breakthrough was obtained by Danchin Eli. However, the regularity and 
the uniqueness of the weak solutions for large data remain largely open for the compress¬ 
ible Navier-Stokes equations, even as in two dimensional space, see Vaigant-Kazhikhov 
m (see also Germain m, and Haspot m, where criteria for regularity or uniqueness 
are proposed). 

The problem becomes even more challenging when the viscosity coefficients depend on 
the density. Indeed, the Navier-Stokes equations (11.311 is highly degenerated at the vacuum 
because the velocity cannot even be defined when the density vanishes. It is very difficult 
to deduce any estimate of the gradient on the velocity field due to the vacuum. This is the 
essential difference from the compressible Navier-Stokes equations with the non-density 
dependent viscosity coefficients. The first tool of handling this difficulty is due to Bresch, 
Desjardins and Lin, see [3], where the authors developed a new mathematical entropy to 
show the structure of the diffusion terms providing some regularity for the density. An 
early version of this entropy can be found in ID for constant viscosity in [351 [36] . The result 
was later extended for the case with an additional quadratic friction term rp|u|u, refer 
to Bresch-Desjardins mm and the recent results by Bresch-Desjardins-Zatorska [5] and 
by Zatorska m- Unfortunately, those bounds are not enough to treat the compressible 
Navier-Stokes equations without additional control on the vacuum, as the introduction of 
capillarity, friction, or cold pressure. 

The primary obstacle to prove the compactness of the solutions to (11.311 is the lack of 
strong convergence for y/~pu in L 2 . We cannot pass to the limit in the term pu (g> u without 
the strong convergence of yfpu in L 2 . This is an other essential difference with the case of 
non-density dependent viscosity. To solve this problem, a new estimate is established in 
Mellet-Vasseur [32] . providing a L°°(0, T; Llog L(fl)) control on p|u| 2 . This new estimate 
provides the weak stability of smooth solutions of (|1.3I) . 
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The classical way to construct global weak solutions of (11.31) would consist in construct¬ 
ing smooth approximation solutions, verifying the priori estimates, including the Bresch- 
Desjardins entropy, and the Mellet-Vasseur inequality. However, those extra estimates 
impose a lot of structure on the approximating system. Up to now, no such approxima¬ 
tion scheme has been discovered. In PI El, Bresch and Desjardins propose a very nice 
construction of approximations, controlling both the usual energy and BD entropy. This 
allows the construction of weak solutions, when additional terms -as drag terms, or cold 
pressure, for instance- are added. Note that their result holds true even in dimension 3. 
However, their construction does not provide the control of the pu in L°°(0, T ; L log L(H)). 

The objective of our current work is to investigate the issue of existence of solutions for 
the compressible Navier-Stokes equations (II.ip with large initial data in 3D. Jungel [24] 

studied the compressible Navier-Stokes equations with the Bohm potential up , and 

obtained the existence of a particular weak solution. Moreover, he deduced an estimate 
of Vp3 in L 4 ((0,T) x fl), which is very useful in this current paper. In [IT] , Gisclon and 
Lacroix-Violet showed the existence of usual weak solutions for the compressible quantum 
Navier-Stokes equations with the addition of a cold pressure. Independently, we proved 
the existence of weak solutions to the compressible quantum Navier-Stokes equations with 
damping terms, see [38] . This result is very similar to m- Actually, it is written in m 
that they can handle in a similar way the case with the drag force. Unfortunately, the 
case with the cold pressure is not suitable for our purpose. 

Building up from the result [38] (a variant of EE), we establish the logarithmic estimate 
for the weak solutions similar to [32]. For this, we first derive a “renormalized” estimate on 
p<^(|u|), for ip nice enough, for solutions of [38] with the additional drag forces. It is showed 
to be independent on the strength of those drag forces, allowing to pass into the limit when 
those forces vanish. Since this estimate cannot be derived from the approximation scheme 
of [38], it has to be carefully derived on weak solutions. After passing into the limit k 
goes to 0, we can recover the logarithmic estimate, taking a suitable function <p. This is 
reminiscent to showing the conservation of the energy for weak solutions to incompressible 
Navier-Stokes equations. This conservation is true for smooth solutions. However, it is 
a long standing open problem, whether Leray-Hopf weak solutions are also conserving 
energy. 

Equation (11.11) can be seen as a particular case of the following Navier-Stokes 
p t + div(pu) = 0 

(pu) t + div(pu <8) u) + VP — div(^(p)Bu) — V(A(p)divu) = 0, J 

where the viscosity coefficients p(p) and A (p) depend on the density, and may vanish on 
the vacuum. When the coefficients verify the following condition: 

A(p) = 2pp(p) - 2 p(p) 

the system still formally verifies the BD estimates. However, the construction of Bresch 
and Desjardins in [3] is more subtle in this case. Up to now, construction of weak solutions 
are known, only verifying a fixed combination of the classical energy and BD entropy (see 
m in the case with additional terms. Those solutions verify the decrease of this so-called 
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K-entropjO, but not the decrease of Energy and BD entropy by themselves. The extension 
of our result, in this context, is considered in [39]. 


The basic energy inequality associated to CEO reads as 

rT r 


m + 


u|“ dx dt < Eq, 


io Jn 


(1.4) 


where 

E(t) = E(p, u )(t) = Qp|u| 2 + p ^ dx, 

and 

E 0 = E(p, u)(0) = Qpo|u 0 | 2 + ^-yP^ dx - 

Remark that those a priori estimates are not enough to show the stability of the solutions 
of (jl.ip . in particular, for the compactness of p 7 . Fortunately, a particular mathematical 
structure was found in mm, which yields the bound of Vpa in T 2 (0, T; L 2 (D)). More 
precisely, we have the following Bresch-Desjardins entropy 

[ f 7 -p|u + Vlnp| 2 H——-^) dx+ [ [ \X7p2\ 2 dxdt 

Jn 7 — 1 / Jo Jn 

+ J J p|Vu - V T u| 2 dxdt < J ^p 0 |u 0 | 2 + |V^/po| 2 + ^ dx. 

Thus, the initial data should be given in such way that 


Po € Po > 0, 

m 0 G T 1 (14), m 0 = 0 


VV^€L 2 (fl), 


if po = 0, 


l m o| 

Po 


g L\n). 


(1.5) 


Remark 1.1. The initial condition V^po G L 2 (Q) is from the Bresch-Desjardins entropy. 


The primary obstacle to prove the compactness of the solutions to (|1.6jl with ro = = 0 

is the lack of strong convergence for ^pu in L 2 . Jungel [2] proved the existence of a 
particular weak solutions with test function p<p, which was used in [4]. The main idea 
of his paper is to rewrite quantum Navier-Stokes equations as a viscous quantum Euler 
system by means of the effective velocity. In [24], he also proved inequality (11,9|i which 
is crucial to get a key lemma in this current paper. Motivated by the works of MEU, 
we proved the existence of weak solutions to (11.61) and the inequality (11.9(1 . see [38]. The 
advantage of ro and r\ terms is that there is a compactness pu (g) u in L 1 and the strong 
convergence of y/pu in L 2 . In particular, we need to recall the following existence result 
in [38]. 


Proposition 1.1. For any n > 0, there exists a global weak solution to the following 
system 

p t + div(pu) = 0, 


{pu) t + div(pu <g> u) + Vp 7 


div(pBu) = —tqu — rip|u| 2 u + KpV( 


A y/p 

Vp 


^Note that k here is not related to the k term in m- 


( 1 . 6 ) 
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with the initial data m and satisfying (USD and —ro f n log_ po dx < oo. In particular, 
we have the energy inequality 

E(t) + f f p\Bu\ 2 dx dt + r 0 f f \u\ 2 dxdt + ri f f p|u | 4 dx dt < Eq, (1.7) 

Jo Jn Jo Jn Jo Jn 

E{t) = E(p, u)(i) = Qp|u | 2 + dx > 


where 


and 


E 0 = E(p, u)(0) = Jn Qpo|u 0 | 2 + ^—[Po + ^l v \/Pol 2 ) dx i 
and the BD-entropy 

j jjp\n + Vln/)| 2 + -?—j + ^\Vyfp\ 2 t 0 log p\ dx+ f f | V p~ 2 1 2 dx dt 
J ri \ O * 1 J J o J 

+ [ f p\ Vu — V T u| 2 dx dt + k f [ p|V 2 log p\ 2 dx dt (1-8) 

Jo Jn Jo Jn 

<2/ (po\n 0 \ 2 + | V-y/pol 2 + ^ 0 1 + — | V-^/po | 2 — ro log_ P(J dx + 2E(j, 


where log_ g = logmin(g,1); the following inequality for any weak solution (p, u) 

K2 \Wp\\L 2 {0,T-,H 2 {n)) + K4 W^P 4 IIl 4 (0,T;L 4 (Q)) < C, 

where C only depends on the initial data. 

Moreover, the weak solution (p, u) has the following properties 

pu € C([0,T\-,Ll eak m, (y/P)t € L 2 ((0,T) x 11); 

If we use (p K ,u K ) to denote the weak solution for k > 0, then 

y/pfu K —>• y/pu strongly in L 2 ((0,T) x 12), as k —)• 0, 
where (p, u) in ED is a weak solution to (HID and (lid for k = 0 . 

Remark 1.2. The energy inequality (11.71) yields the following estimates 

llv / P u llL°°( 0 ,T;L 2 (n)) < E 0 < OO, 
l|p||Z/ oo (0,T;LT'(n)) < Eq < OO, 

\\VKVy/p\\L°°(0,T-,L 2 (n)) < E 0 < OO, 

llv / /® U IU 2 ( 0 ,T;L 2 (n)) < Eo < 00 , 
ll\/^ 0 u IU 2 ( 0 ,T;L 2 (n)) < Eo < 00 , 

H^lP u llL 4 (0,T;L 4 (r2)) < ^0 < OO. 

The BD entropy (11.81) yields the following bounds on the density p: 

ll^ 7 VPllL°°( 0 ,T;L 2 (n)) <C<00, 
|| v ^cpV 2 logp|| L 2 ( 0 )T;i 2 (n)) < C < oo, 

7 

l|V/5 2 IU 2 ( 0 ,T;L 2 (n)) < C < OO, 


(1.9) 


( 1 . 10 ) 

( 1 . 11 ) 


( 1 . 12 ) 


(1.13) 

(1.14) 

(1.15) 
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and 

ll\/P Vu llL 2 ( 0 ,T;L 2 (n)) < (po|u 0 | 2 + + |V^/Po | 2 - r 0 log. Po^j dx + 2 Eq < oo, 

n.i(>) 

where C is bounded by the initial data, uniformly on ro, r i and k. 

In fact, (|1.13D yields 

^p€L°°(0,T-,L 6 (Q)), (1.17) 

in three dimensional space. 

Remark 1.3. Inequality (11.91) is a consequence of the bound on (11.141) . This was used 
already in [23] . The estimate for the full system (11.61) is proved in 


Remark 1.4. The existence result of [2J contained the case with k = 0, which can be 
obtained as the limit when k > 0 goes to 0 in (11.61) . by standard compactness analysis. 

Remark 1.5. The weak formulation reads as 


pu ■ iJj dx\l = Q — / / puv/y dx dt — / / pu <g) u : S7if dx dt 

Jo Jn Jo Jn 

[ I p 7 divi/' dx dt — f I pDu : V if’dx dt 

Jo Jn Jo Jn 

T 


r T 

iq Jn 
= ~ro 


(1.18) 


uijj dx dt — n 


io Jn 


/o Jn 


p|u| 2 ui/’ dx dt — 2 k 


/o Jn 


A^JpV y/ptp dx dt 


— K 


[ [ A VP\/pd 
Jo Jn 


i v'lp dx dt. 


>o Jn 
for any test function i/>. 

Our first main result reads as follows: 


Theorem 1.1. For any 5 € (0,2), there exists a constant C depending only on 5, such 
that the following holds true. There exists a weak solution (p, u) to (USD with k = 0 
verifying all the properties of Proposition 17771 and satisfying the following Mellet-Vasseur 
type inequality for every T > 0, and almost every t < T: 


/ p(f, x)(l + |u(t, x)| 2 ) ln(l + |u(t, x)\ 2 ) dx 

Jn 

< J p 0 (l + |u 0 | 2 ) ln(l + |u 0 1 2 ) dx + 8 J ^p 0 |u 0 | 2 + ^y+ |V % /p^| 2 -r 0 log_p 0 ^ dx 

+ 16E 0 + C [ f [ (p 27-1- ^)^ J ( [ p(2 + ln(l + |u| 2 ))« dx ) dt, 

Jo \Jn J \Jn J 


where 7 > 1 in two dimensional space and 1 < 7 < 3 in three dimensional space. 


Remark 1.6. The right hand side of the above inequality can be bounded by the initial 
data. In particular, it does not depend on ro and n. This theorem will yield the strong 
convergence of in space T 2 (0, T ; Q) when ro, ri converge to 0. It will be the key tool 
of obtaining the existence of weak solutions, in [32] . 
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We define the weak solution (p, u) to the initial value problem (11.11) in the following 
sense: for any t E [0, T], 

• ( 11 . 21 ) holds in T>'(Q), 

• (11.41) holds for almost every t € [0, T], 

• m holds in T)' (({), T) x 0)) and the following is satisfied 

P> 0, P EL~([0,T];L^)), 

p(l + |u| 2 )ln(l + |u| 2 ) E L°°(0,T; L 1 (J7)), 

Vpi E L 2 (0, T; L 2 (fl)), V,/p E L°°(0, T; L 2 (0)), 

^/pu E L°°(0, T; L 2 (fl)), /pVu E L 2 (0, T; L 2 (ft)). 

Remark 1.7. The regularity Vy/pE L°°(0, T; L 2 (fl)) and Vpz E L 2 (0, T; T 2 (fi)) are from 
the Bresch-Desjardins entropy. 

As a sequence of Theorem 11.11 our second main result reads as follows: 

Theorem 1.2. Let (po, mo) satisfy (11.51) and 

/ po(l + |u 0 1 2 ) ln(l + |u 0 | 2 ) dx < oo. 

Jn 

Then, for 7 > 1 in two dimensional space and 1 < 7 < 3 in three dimensional space, and 
any T > 0, there exists a weak solution of (pTp-CLZD on ( 0 , T). 


We cannot obtained directly the estimate of Theorem 11.11 from (11.61) with k = 0, because 
we do not have enough regularity on the solutions. But, the estimate is not true for the 
solutions of (ll.fip for k > 0. The idea is to obtain a control on 

/ p(t,x)<p n ((/)(p)u(t,x))dx 

Jn 

at the level k > 0, for a ip n , suitable bounded approximation of (1 + |u| 2 )ln(l + |u| 2 ), 
and a suitable cut-off function <f of p, controlling both the large and small p. The first 
step (see section 2) consists in showing that we can control (uniformly with respect to k) 
this quantity, for any weak solutions of (11.61) with k > 0. This has to be done in several 
steps, taking into account the minimal regularity of the solutions, the weak control of the 
solutions close to the vacuum, and the extra capillarity higher order terms. In the limit 
k goes to zero, the cut-off function f> has to converge to one in a special rate associated 
to k (see section 3 and 4). This provides, for any weak limit of (11.61) obtained by limit n 
converges to 0, a (uniform in n, tq, and 77) bound, to: 


p(t, x)ip n ((j)(p)u(t, x)) dx. 


Jft 

Note that the bound is not uniform in n, for k fixed. However, it becomes uniform in n at 
the limit n converges to 0. In section 5, we pass into the limit n goes to infinity, obtaining 
a uniform bound with respect to rg and rq of 


p(f, x)(l + |u (t, x)| 2 ) ln(l + |u(f, x)\ 2 ) dx. 


Section 6 is devoted to the limit r 1 and ro converges to 0. The uniform estimate above 
provides the strong convergence of ^/pu needed to obtain the existence of global weak 
solutions to a with large initial data. 
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2. Approximation of the Mellet-Vasseur type inequality 

In this section, we construct an approximation of the Mellet-Vasseur type inequality for 
any weak solutions at the following level of approximation system 

Pt + div(pu) = 0, 

o A ,fp (2.1) 

(pu) t + div(pu (g) u) + Vp 7 — div(pBu) = — rou — rrpluru + np\7( - ), 

s/P 

with the initial data (USD, verifying in addition that po > for mo > 0 and ^/pqvlq € 
L°°(0). This restriction on the initial data will be useful later to get the strong convergence 
of s/pu when t converges to 0. This restriction will be cancel at the very end, (see section 
6 )- 

In the same line of Bresch-Desjardins El SUM], we constructed the weak solutions to 
the system (11.611 for any k > 0 by the natural energy estimates and the Bresch-Desjardins 
entropy, see [38] • The term rou turns out to be essential to show the strong convergence of 
s/pu in L 2 (0, T; L 2 (fl)). Unfortunately, it is not enough to ensure the strong convergence 
of yfpu in T 2 (0, T; L 2 (D)) when r$ and r\ vanish. 

We define two C°°, nonnegative cut-off functions f> m and fix as follows. 

<t>m(p) = 1 for any p > —, f m {p) = 0 for any p < (2.2) 

m 2m 

where m > 0 is any real number, and < 2m; 

and 4>k{p) £ C 00 ^) is a nonnegative function such that 

4>k{p ) = 1 fo r any p < K, 4>k(p) = 0 for any p > 2 1\, (2-3) 

where K > 0 is any real number, and 

We define v = <f>(p) u, and <j){p ) = (f rn {p) ( t ) K{p)- The following Lemma will be very useful 
to construct the approximation of the Mellet-Vasseur type inequality. The structure of 
the k quantum term in |23] is essential to get this lemma in 3D. It seems not possible to 
get it from the Korteweg term of BD [2] in 3D. 

Lemma 2.1. For any fixed, k > 0, we have 

II^ v IIl 2 (0,T;L 2 (T2)) < Cj 

where the constant C depend on k > 0, r\, K and m; and 

pt £ L 4 (0, T ; L 6 / 5 (D)) + L 2 (0, T ; L 3 / 2 (0)) uniformly in n. 

Proof. By (11.91) . we have 

ll^ 7 / 94 IIl 4 (0,T;L 4 (D)) A C. 

For v, we have 

Vv = V(^(p)u) = (0'(p)Vp) u + 0(p)V u, 

and hence 

||(^'(p)Vp)u + </ , (p)Vu||^ 2 (0j 2 ’ ; / / 2 (f 2 )) 

< Cj|p 4 uVp 4 ||l 2 ( 0 , T ; L 2 (D)) + C < II\/P^ U IIl 2 ( 0 ,T;I, 2 (D)) 

< C'||p 4 u|| L 4( 0iT;i 4(Q)) ||Vp4 11L 4 ( 0 ,T ; L 4 (Q)) + C|l\/P^ U llL 2 (0,T;L 2 (n))) 
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where we used the definition of the function fi(p). Indeed, there exists C > 0 such that 

<Kp) 


\4>'{p)^/p\ + 


VP 


< c 


for any p > 0 . 

For pt, we have 

p t = — V /9 • u — pdivu 

= — 2 Vy fp- piupi — -y/p-y/jodivu = Si + S' 2 . 
Thanks to (| 1 . 12 |) . (11.131) and (11,171) . we have 


Si € L 4 ( 0 ,T;L r (H)) for 1 < r < 


6 


By (|1 . 12|) and (11.171) . we have 


Thus, we have 


S 2 € L 2 ( 0 ,T;L s (H)) for 1 < s < 


p t E T 4 (0, T; L r (fi)) + L 2 (0, T; L s (fi)). 


We introduce a new C' 1 (M 3 ), nonnegative cut-off function p n : 

<p n (u) = (Pn(|u| 2 ), 

where <p n is given on M + by 

if 0 < y < n, 
if n < y < C n 


finiv) < 


_ l 
1 + 2 / 

__I_ 

1 + 2 / 

= 0 if y > C n , 


with (p' n ( 0) = 0, (p n (0) = 0, and C n = e(l + n) 2 — 1. 

Here we gather the properties of the function (p n in the following Lemma: 

Lemma 2.2. Let ip n and <p n be defined as above. Then they verify 

• a. For any u E M 3 , we have 

Tni u ) = 2 ( 2< /"(l u | 2 ) u ® u + I( £n(l u | 2 )) , 

where I is 3 x 3 identity matrix. 

• b. /or any n > 0 and any y > 0. 


c. 


^n(y)< 


□ 


(2.4) 


(2.5) 


( 2 . 6 ) 


(2.7) 


= 1 + ln(l +y) if0<y<n , 

= 0 ify>C n , 

> 0, and <l + ln(l + y) if n < y < C n . 

In one word, 0 < (p' n < 1 + ln(l + y) for any y > 0, and it is compactly supported, 
d. For any given n > 0, we have 

\<Pn( u )\ < 6 + 21n(l + n) (2.8) 

for any u £ tl 
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• e. 


<Pn(y) 


(1 + y) ln(l + y) if 0 < y < n, 

< 2(1 + ln(l + n))y — (1 + y)ln(l + y) + 2(ln(l + n) — n), if n<y < C ni 
e(l + n) 2 - 2n - 2 if y > C n , 


(2.9) 

<p n (y) is a nondecreasing function with respect to y for any fixed n, and it is a 
nondecreasing function with respect to n for any fixed y. 


<Pn{y) -1 (1 + y) ln(l + y) a.e. 


( 2 . 10 ) 


as n —>• oo. 


Proof. We prove each statement one by one as follows: 

• a. Thanks to (12.41) . we have ) = 2</3( l (|u| 2 )u, and 


+n( u ) = 2 ( 2 < ^n(l u | 2 ) u ® U + I^(|U| 2 )) , 


where I is 3 x 3 identity matrix. 

• b. The statement of b. follows directly from (12.51) . 

• c. Integrating (12.5p with initial data ^(0) = 0, we have 


p'niy) 


1 + ln(l + y) if 0 < y < n, 

< 1 + 21n(l + n) - ln(l +y), if n < y < C n 
^0 if y > C n , 


Since 

1 + 2 ln(l + n) — ln(l + C n ) = 0, 
thus, for any y > 0 , we have 


p'niy) > 0 - 


For any n < y < C n , we have 


( 2 . 11 ) 


1 + 21 n(l + n) - ln(l + y) < 1 + 21 n(l + y) — ln(l + y) = 1 + ln(l + y). 

Thus, for any y > 0, we have 

<PM < 1 + ln(l + y). 

• d. By a.-c., 

|+n( u )l < 4 l+nll u ! 2 + 2 |+nl < ["! , 2 + 2(1 + ln(l + n)) < 6 + 21 n(l + n). 

• f. Integrating (12.111) with initial data ip n ( 0) = 0, it gives (12.91) . Moreover, thanks 
to c., <Pn(y) is an increasing function with respect to y for any fixed n. We have 
also that <p n (y) is a nondecreasing function with respect to n for any fixed y. 

□ 


The first step of constructing the approximation of the Mellet-Vasseur type inequality 
is the following lemma: 
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Lemma 2.3. For any weak solutions to (EH) constructed in Proposition \1.1[ and any 
ijj(t) € £>(— l,+oo), we have 

- [ [ iptPVniy) dx dt + f f ip(t)ip' n (y)F dx dt + f f ip(t)S : V(<p' n (v)) dxdt 

Jo Jn Jo Jn Jo Jn 

= / p 0 <p n (v 0 )il}(0)dx, 

Jn 


where 


S = p0(p)(Bu + k ——I), and 

Vp 


F = p 2 U(f>' (p)divu + 2 p 2 Vp 2 0(p) + pV0(p) Du + r 0 uc/)(p) 

+ np|u| 2 u(/Hp) + Ky/pV(/)(p)A^p + 2 kc/)(p)'V y/pAy/p, 

where I is an identical matrix. 


( 2 . 12 ) 


(2.13) 


In this proof, k, m and I\ are fixed. So the dependence of the constants appearing in 
this proof will not be specified. 

Multiplying 0(p) on both sides of the second equation of (|2.1I) . we have 
(pv) t - pu<j/(p)p t + div(pu <g) v) — pu ® uV 0 (p) + 2 pa Vp 2 0 (p) 

- div(0(p)pBu) + pV^(p)Du + r 0 u</>(p) + np|u| 2 U(/>(p) - «V(- y /p0(p)A A /p) 

+ KyfpS7(j){p)Ayfp + 2k4>(p)V y/pAyfp = 0. 

Remark 2.1. Both V^/p and pf are functions, so the above equality are justified by regu¬ 
larizing p and passing into the limit. 

We can rewrite the above equation as follows 

(p v )i + div(pu ( 8 ) v) — divS + F = 0, (2.14) 

where S and F are as in (12.131) . and here we used 

pucj)'(p)p t + pu (8) u 4>{p)Vp = p\i4>{p)(p t + Vp • u) 

= — p 2 U(//(p)divu. 

We should remark that, thanks to (11.91) . (|1.12D - (I1.16I) . 

since y/pcj>(p) and p<J>(p) bounded. Those bounds depend on K and k. 


We first introducing a test function ij;(t) € 20(0, +oo). Essentially this function vanishes 
for t close t = 0. We will later extend the result for ip(t) € 2)(—1, +oo). We define a new 
function $ = where f(t,x) = f * ?/&(£, x), k is a small enough number. Note 

that, since ij){t) is compactly supported in (0,oo). 4? is well defined on (0,oo) for k small 
enough. We use it to test (|2.14l) to have 

f [ ii’(t)ip' n (y){(pv) t + div(pu<8> v) — divS + F] dxdt = 0, 

Jo Jn 
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which in turn gives us 

rT 


f f 'ip(t)(p' n (v)[(pv) t + div(pu <S> v) — divS + F] dxdt = 0. 

Jo Jn 

The first term in (12.151) can be calculated as follows 
/ / i’(i)(p , n(v)(pv) t dxdt 

Jo Jn 

= [ [ i’(,t)ip' n (y)(pv)tdxdt+ [ f il>(t)(f/ n (v)[{pv)t-(pv)t]dxdt 

Jo Jn Jo Jn 


(2.15) 


+ R\ 


where 


/ / '<P{t)<p , n ('v)(ptv + pv t ) dx dt 

Jo Jn 

/ / i/j(t)p t (p' n {V)Vdxdt+ / / ^(t)pp n (w) t dxdt +1 

Jo Jn Jo Jn 

[ [ *l>(t)<Pn(v)[(pv)t - (pv)t]dxdt. 

Jo Jn 


(2.16) 


Ri = 


Thanks to the hrst equation in (12.11) , we can rewrite the second term in (12.151) as follows 

rT 


[ [ 'll)(t)<p' n (v)div(pu®v) dxdt 

Jo Jn 

iit)p t ip n (v)dxdt - / / il>(t)p t ip' n (y)v + Rq, 

i Jo Jn 


f 

/o Jn 


(2.17) 


and 


R2 = 


to Jn 


■0(t)</? n (v)[div(/9U (g> v) — div(pu ® v)]. 


By (|2.15j) - (|2.17j) . we have 

rT 


[ [ V’( i )(P 7 ? n(v))t dxdt + R\ + R 2 - f f V , (i)< 7 , n( v )divS dxdt 

Jo Jn Jo Jn 

rT 


+ 


[ [ ^(t)ip'n{-v) F = 0 . 

Jo Jn 


(2.18) 


Notice that v converges to v almost everywhere and 

pip n (y)il) t -A pip n (y)i> t in L\(0,T) x 0). 

So, up to a subsequence, we have 

/ / [piPn(y : ))'4>t dx dt —>• / / {pip n {'v))ipt dx dt as k —» 0. (2-19) 

jo Jn Jo Jn 

Since ^(v) converges to ip' n {y) almost everywhere, and is uniformly bounded in L°°(0, T; 12), 
we have 


Noticing that 


[ [ J>{tW n (y)F -A f [ l/)(t)(fJ n (v)F 

Jo Jn Jo Jn 


as /c —»• 0. 


( 2 . 20 ) 


Vv g L 2 (0,T-L 2 (Q)), 
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we have 

Vv-> Vv strongly in T 2 (0, T; L 2 (fl)). 

Since S converges to S strongly in T 2 (0, T; T 2 (17)), and ’p'f (v) converges to <p”(v) almost 
everywhere and uniformly bounded in L°°((0, T) x 17), we get 


f f -i/>(t)<// n (v)divS dx dt = — 
to ./n 

which converges to 

~[ [ if{t)S :V(ip' n {v)) dxdt. 

Jo Jn 

To handle R\ and R 2 , we use the following lemma due to Lions, see 

Lemma 2.4. Let f E PL 1,p (M iV ), g € L 9 (M Ar ) with 1 < p, q < 00 , and ^ | < 1. Then, 

we have 


[ if(t)S : V(</?' n ( v )) dxdt, (2.21) 
Jn 

( 2 . 22 ) 


l|div(/ 5 ) *w e -div(/( 5 *w e ))|| L r (R iv ) < C'||/|| w i, P(R jv ) || 5 r|| i9(K v ) 
for some C > 0 independent of e, f and g, r is determined by £ = | + |. /n addition, 

div(fg) *w e - div(f(g * w e )) 40 in TfK^) 


1 _ 1 1 1 

r 3 
r t /n3>iV\ 


as e -4 0 if r < 00 . 


This lemma includes the following statement. 

Lemma 2.5. Let /* € T p (0,T), g E L 9 (0,T) with 1 < p, q < 00 , and | | < 1. Then, 

ice have 

IK/flOt * We — {f (9 * We))t\\L r (0,T) < C11 /i 11 Lr> (0,T) 11511M (0,T) 

/or some (7 > 0 independent of e, f and g, r is determined by £ = | + |. In addition, 

( fg)t *w e - (f(g * w e )) t -4 0 in T r (0, T) 

as e ^ 0 if r < 00 . 


With Lemma [2.41 and Lemma [2.5 1 in hand, we are ready to handle the terms and i? 2 . 
For k > 0, by Lemma 12.II and Poincare inequality, we have v E T 2 (0, T; L 6 (!2)). We also 
have, by Lemma 12. 11 

pt € T 4 (0,T;L 6 / 5 (L!)) + T 2 (0,T;L 3 / 2 (L!)). 


Thus, applying Lemma YI7A 

,-T 


\Ri\< f [ il>(t)<ft n (v)[(pv)t - (pv)i] 
Jo Jn 

<C(if) [ [ (f/ n (v)[(pv) t - (pv) t ] 

JO ill 

Similarly, applying Lemma 12.41 we conclude 


dx dt 

dx dt -4 0 as k -4 0. 


(2.23) 


R 2 ^ 0 as k —^ 0. 


(2.24) 
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By (I2.19M2.24D . we have 


n 


iptP<Pn(v) dxdt + 


lo Jn 


ip(t)ip' n (v)F dx dt 


+ 


/ f ^{t)S : V{ip' n (v)) dx dt = 0, 

Jo Jn 


(2.25) 


for any test function -0 € 0(0, oo). 

Now, we need to consider the test function ip(t) £ D(—l,oo). For this, we need the 
continuity of p(t) and (y / pu)(t) in the strong topology at t = 0. 

In fact, thanks to Proposition ll.il we have 


(Vp)t € L 2 ( 0 , T; L 2 m, yfp £ L 2 ( 0 , T; H 2 (Q)). 


This gives us 

yfpe C([0,T];L 2 (11)) and Vy/p £ C(0, T; L 2 (H)), 


thanks to Theorem 3 on page 287, see [12] . Similarly, we have 


p£C([0,T];L 2 (F!)) 


(2.26) 


due to 

1 3 

II^ 7 P||l 2 (0,T;L 2 (O)) < C , ||V/9 :5 || i 4( 0jT . L 4(f2))||p4 || i4 ( 0)T . L 4(Q)). 
Meanwhile, we have 


sjp £ L°°(0, T; L P (Q)) for any 1 < p < 6, 


and hence 


\fp £ C([0, T]; L P (Q)) for any 1 < p < 6. 
On the other hand, we see 


(2.27) 


ess 


limsup / \y/pu - VPouol 2 dx 
t —>0 J £1 


7 

Po 


< ess limsup / (-p|u| 2 H -- + n\Vy/p\ 2 ) dx - / (-poK| S-) + k|Vvp 0 

t->o V ./o 2 Y — 1 ./o 2 Y — 1 


' dx 


7 


7‘ 




+ ess limsup ( 2 / ^/^u 0 (Vpou 0 - VP U ) dx + ( 


7 

Po 


7 — 1 7 — 1 


— k ess lim sup | V yfp Q — V yfp\ dx + 2k ess lim sup / V y/p 0 ■ (V y/p Q — V y/p) dx. 

t—>0 t—>0 J £1 


(2.28) 


We have 


ess lim sup 

t—>0 , 


'n 


v \/p 0 • ( v \/Po - VVp) dx = 0. 


(2.29) 
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So, using (11.71) . (12.271) and the convexity of p i-a p 7 , we have 
esslimsup / IVP u -\/A) u ol 2 dx 

t —>-0 J £7 


< 2ess 


limsup / VPo u o(\/A)U 0 - yfp\i)dx 
t—J r2 




= 2ess limsup ( / VPo u o(V~Po u o - ^puc/) m (p)) dx + / v /p^u 0 (l - (j) m (p))^pudx 
t — \Jq 

= B\ + i?2- 

By Proposition ll.il we have 


pueCaO.T];^^)). 


We consider B± as follows 
B 1 = 


2esslimsup( [ v / Po u o( ^t^ (AW o ~ pu)) da; ~ [ y/PoPo\ u o\ 2 ( 
t->o \Jn VP Jn 

then we have B\ = 0, where we used (12.271) and (12.301) . 


fimip) $m(Po)' 

Vp Vpo ' 


(2.30) 


dx 


Since m > mo, and po > we have 

1-^21 < llv / PO U ollL 00 ( 0 ,T ; a)llv / P U llL 00 ( 0 ,T;I, 2 (Q)) ess li msu P P — 4>m (p) II L 2 ( 0 ,T;O) = 0. 

i —>0 


Thus, we have 


ess 


which gives us 

By (I2.26p and (12.311) . we get 

lim — 

t— >o r 


/o in 

Considering (|2.25p for the test function, 


limsup / |Ppu - VPou 0 | 2 da; = 0, 
t —>-0 J r2 


VpugC([0,T];L 2 (P!)). 


pip n (v) dx dt = / po^n(vo) dx. 


(2.31) 



VtV) = ip(t) for t > T, lj) T (t) = V{ T )— f° r t < T, 


we get 

rT 


- / VtPPniy) dx dt + / / 1 pr(t)ip'ni V )Fdxdt 

J r J Q J 0 J Q 

+ [ [ Vv(i)S : V(y 4 (v)) dxdt — ^ 

Jo Jn 

Passing into the limit as t —>• 0, this gives us 

- [ [ VtP'Pn (v) dx dt T [ [ V(t)<p' n {v)F dxdt 

Jo Jn Jo Jn 



T Jo Jn 


pip n (v) dx dt. 


+ 


f ( ip(t)S : V(ip'n(v)) dx dt = j poi>(0)ip n (y 0 ) dx. 

Jo Jn Jn 


(2.32) 
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3. Recover the limits as m — >• oo 

In this section, we want to recover the limits in (12.121) as rn oo. Here, we should 
remark that (p, u) is any fixed weak solution to (12.11) verifying Proposition 11.11 with k > 0. 
For any fixed weak solution (p, u), we have 

4>m{p ) —> 1 almost everywhere for (t,x), 

and it is uniform bounded in L°°(0,T;R), we also have 

r 0 feWueI 2 (0,T;L 2 (O)), 

and thus 

v m = (j)m4>K u —> c/^u almost everywhere for (t, x ) 
as m —> oo. By the Dominated Convergence Theorem, we have 

v m ->• 0a: u in L 2 (0,T;L 2 (R)) 
as m —> oo, and hence, we have 

^n(v m ) ->■ ^ n (0A'u) in L p ((0,T) x R) 


for any 1 < p < oo. 

Meanwhile, for any fixed p, we have 

(j>' m {p) 0 almost everywhere for (t,x) 

as m —>• oo. We calculate \<f' m (p)\ < 2m as ^7 < p < 7 ^, and otherwise, (f>' m {p) = 0, thus, 
we have 

IM™(P)I < 1 for a11 P- 


We can find that 
r T 


/ / il>'{t)(p<p n (v Tn ))dxdt -A / / if 1 {t){pip n {<i>K{p)v)) dx dt 

J 0 1!) JO J!) 


and 


/ POPn(Vmo) ^ / PO < Pn(/ ) Ar(Po)uo) 
JSI J[J 


asm-> 00 . 

To pass into the limits in (12.321) as m -> 00 , we rely on the following Lemma: 

Lemma 3.1. If 


||||L oo (o,T';fi) < C, a m —^ a a.e. /or (t, x) and in L p ((0,T) x R) /or any 1 < p < 00 , 
/ E L 1 ((0, T) x R), f/ien we hare 


/ / <t>m(p)ci m f dx dt -A / / afdxdt 

Jo Jn Jo Jn 


as m -A 00, 


and 
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Proof. We have 

\fm{p)a m f ~ a/| < \(p m {p)f ~ f ||a m | + |a m / - a/| = I\ + h- 
For Ip. 4> m {p)f —> / a.e. for (t,x) and 

\<t>m{p)f ~ f\< 2|/| a.e. for (t,x), 
by Lebesgue’s Dominated Convergence Theorem, we conclude 


f f \<t>m{p)f ~ f\dxdt -i 0 
JO Jo 


asm-loo, which in turn yields 


|0m(p)a m / - a m /| dxdf 


< ||a m ||z,°°(0,T;O) / / I $m{p)f /1 dx dt 

JO JO 


as m —l oo. Following the same line, we have 


[ [ | a m f ~ a/| dx dt-> 0 

Jo Jo 


as m —> oo. Thus we have 


( f > m{p)a rn f dx dt 


af dx dt 


as m —l oo. 


We now consider / Q T f n \p(f>' m (p)a m f\ dx dt. Notice that |p</>( n (/9)| < (7, and p<p' m {p) con¬ 
verges to 0 almost everywhere, so \pc/)' m (p)a m f\ < C\f\ , and by Lebesgue’s Dominated 
Convergence Theorem, 

/ / lMn(p)a m /l dxdt^-0 

Jo Jo 

as m —l oo. □ 


Calculating 


/(/S m : V (ip' n (v m )) dx dt 


where 


'0 Jo 

rT r 


^{t) 

^m^Pnfarri ) (V</> m / A -U + 0 m V/ A 'U + / m / A -Vu) dx dt (3.1) 

4>m{p)a m ifidxdt+ [ [ p4> rr fp)a na f- 1 dxdt, 


Clm 1 = 4>m{p)Pn{y m ), 

TTTl,, l ... ^y/~P tt\ /,,Y7JL. 


/ 1 = il>{t)p<j> K (p) + A—— Ij (uVfe(p) + ^jf(p)Vu), 


a m 2 = 9 ? n( v m)<( , m(p)0A-(/O)u = <^"(v m )v m , 
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h = i/}(t)(j>K(p) ^Du + K — V P = 2'4(t)4 K {p)(KAy/pVy/p + y/pDuVy/p). 
So applying Lemma 13.11 to (13.11) , we have 

[ [ ip(t)S m : V (<Pn{v m )) dx dt -A [ [ ip(t )S : V(^' n (^(p)u)) dx dt 

Jo Jn Jo Jn 

as m -A oo, where S = ^(p)p(Du + k^£-T). 

Letting F m = F m i + F m2 , where 

F m i = p 2 U(//(/9)divu + pS74{p) Du + KyfpSJ <f>(p) Ay/p 

a Vp, 


where 


and 


where 


= P {4m(p)4K(p) + </>m(p)</4(p)) (pudivu + Vp • Du + KVp— 
(/>A'(p)(p u divu + Vp • Du + kV p——) € L 1 ((0, T ) x 0), 

Vp 

p4'kV ')(pudivu + Vp • Du + «Vp——) € L 1 ((0, V) x 0), 

Vp 


F m 2 — 4 m (p)4k(p)Vp^ 2 Vp 2 + r 0 u + np|u| 2 u + 2FS7yfp/S.yfp), 


where 


4 k (p) 2 Vp 2 + r 0 u + rip|u| 2 u + 2 kVvpAVP^ € L 1 ((0,T) x 0). 

Applying Lemma I3.1L we have 

[ [ 4>(t)ip' n (v m )Fm dxdt -A [ ( 4(t)ip' n (4K(p)u-)F dxdt, 

Jo Jn Jo Jn 

F = p 2 u4' K (p)dwu + 2 p^Vp^4 K (p) + pV^A:(p)Du + r o u0A'(p) 

+ rip|u| 2 u 0 A'(p) + kVp'V 4 k(p)&Vp + 2 k 4 k(p)V Vp A Vp- 

Letting m —> oo in (12.321) , we have 

- [ [ VV){P l Pn(4K{p)V)dxdt+ [ [ 4(t)Vn(4K(p)V F dxdt 

Jo Jn Jo Jn 

+ / / i>{t)S : V(<p' n {4K(p)u))dxdt = / ■4(0)p 0 (p n (4 K (p 0 )u 0 ) dx, 

Jo Jn Jn 

which in turn gives us the following lemma: 

Lemma 3.2. For any weak solutions to m verifying in Proposition \l.ll we have 

- [ [ '4 , (t){pPn(4K(p)V)dxdt+ [ [ -4(t)ip'n(4K(p)u)F dxdt 

Jo Jn Jo Jn 

+ / 4(t)S :V(Vni4K(p)V)dxdt = / 4(0)p 0 (p n (4 K {po)uo) dx, 

Jo Jn Jn 
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where S = ^(p)p(Du + and 

F = p 2 u(/>' K (p)dwu + 2p2 Vp%4>k(p) + pV0k(p) Ou + r 0 u/x(p) 

+ np\u\ 2 U(f)K(p) + K^p\7(/)k(p)A^P + 2K<j> K {p)V yfp&y/p- 
where I is an identical matrix. 

4. Recover the limits as k -a 0 and K -a oo. 

The objective of this section is to recover the limits in (13.21) as k -a 0 and K —> oo. In 

_ 3 

this section, we assume that K = k 4 , thus /i -A oo when k -a 0. First, we address the 
following lemma. 

Lemma 4.1. Let k -a 0 and K -A oo, and denote v K = <I>k(Pk) u«, 
pZ is bounded in L r ((0,T) x D) /or any 1 < r < 2 
5 

and any 1 < r < - in 3D. 

For any g E C' 1 (M + ) with g bounded, and 0 < a < oo in 2D, 0 < a 
PkP(M 2 ) p“p(|u| 2 ); strongly in L l ((0,T) x 
In particular, we have, for any fixed n, 

Purify k) PPn(u) strongly in L 1 ((0,T) x II), (4.2) 

and 

p 27_L (l + ( Pn(M 2 )) -> P 27 "^ 1 + <Pn(l u | 2 )) strongly in L 3 (( 0, T) x 14), (4.3) 

for 1 < 7 < 3. 

Proof. In 2D, we deduce that 

pZ € L°°(0, T : L 1 (D)) D L^/O, T; L P (D)) for any 1 < p < oo. 

Thus, pZ is bounded in L r (( 0 , T) x D) for 1 < r < 2 . 

In 3D, we deduce that 

P 7 € i°°(0, T : L\Q)) n ^(0, T; L 3 (D)). 

Applying Holder inequality, we have 

2 3 

Thus, pZ is bounded in L3((0,T) x D). 

We have that (p K )t is uniformly bounded in 

L 4 (0, T; L 6 / 5 (D)) + L 2 (0, T; L 3 / 2 (D)), 
thanks to Lemma ED and also we have 

||Vp /t || L cx.(o j T;L 3 / 2 (n)) ^ 

Applying Aubin-Lions Lemma, one obtains 

p K -A p strongly in L p { 0, T ; L 3 / 2 (I4)) for p < oo. 


u>e have 
in 2D, 

(4.1) 

< ^ in 3D, we have 
fi). 
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When k —^ 0, we have u K —> yfpu strongly in L 2 (0, T; L 2 (fl)) in Proposition 11.11 ( 
also see [38]). Thus, up to a subsequence, for almost every (t,x) such that p(t,x) ^ 0, we 
have 

u K (t,x) = ——-)• u(t,x), 

v P K 

and 

v K ->• u(f,x), 

as k -A 0. For almost every (t, x) such that p(t, x ) = 0, 

|p“ 5 (I v k| 2 )| < Cp*{t,x) -A 0 = p°5(|u| 2 ) (4.4) 

as k —^ 0. 

Hence, pfg (|v re | 2 ) converges to p Q g , (|u| 2 ) almost everywhere. Since g is bounded and (14.111 . 
p°g(\v K \ 2 ) is uniformly bounded in L r ((0,T) x Q) for some r > 1. Hence, 

P°9( lw| 2 ) -> P“ff(|u| 2 ) in T 1 ((0,T) x H). 

By the uniqueness of the limit, the convergence holds for the whole sequence. 

Applying this result with a = 1 and g( |v re | 2 ) = ip n (y K ), we deduce (14.21) . 

Since 7 > 1 in 2D, we can take a = 2q — 1 < 2q; and take 7 < 3 in 3D, we have 
27 — 1 < Thus we use the above result with a = 27 — 1 and g{ |v K | 2 ) = 1 + <Pn(|v K | 2 ) 
to obtain (14.31) . 

□ 


With the lemma in hand, we are ready to recover the limits in (13.21) as n 
K —>• 00. We have the following lemma. 

_ 3 

Lemma 4.2. Let K = k 4 , and k -a 0, for any if > 0 and ip < 0, ioe have 

- / ip'(t)pcp n (u)dxdt 

Jo Jft 

< 8||V>||l°= (Yl u ol 2 + + l v \/Pol 2 - r olog_ + 2 f; c 

+ CdIV’IU 00 ) [ [ ( 1 + ( Pn(\ u \ 2 ))p 2l ~ 1 dxdt + 'ip( 0) [ p 0 <p n (u 0 )dx 

Jo Jft Jft 


0 and 


(4.5) 


Proof. Here, we use (p K , u K ) to denote the weak solutions to (12.11) verifying Proposition 
11.11 with k > 0. 

By Lemma l4.ll we can handle the first term in (13.21) . that is, 


and 


[ [ ij/(t)(p K (p n (y K ))dxdt -A [ [ p/ (t)(pip n (u)) dxdt 

Jo Jft Jo Jft 

ip(0) / Po^'„( v «,o) dx -A ^(0) / pQ(p' n ( u o) dx 


as K —>■ 0 and K = k 4 -A 00. 


(4.6) 


(4.7) 
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V K ) ■'Vp 1 K (j) K (p K )dxdt = 


[ [ 

Jo Jn 

~ / , <P(t)pZ<l> K (p K ,)<Pn : 'Vv K dxdt- / / il>(t)pl<p' n (v K ) ■ V<t> K (p K ) dxdt 

Jo Jn Jo Jn 


= Pl+P2- 

We can control P\ as follows 


\P“2 1 — 


L°° 


>0 Jn 


IPlWnW^^Kip^ldx dt 


<C(n, ||V’||l°°)« *\\(I>'kVp^\\l°°\\p, 


7+3, 


< -J=C(n, HVIU~)« i ll^ +i ll i i (0iT .. i # (SI)) 


L3(0,T;L3(Q)) 
1 


N 4 || yP k ||l 4 (0,T;L 4 (0)) 
N 4 || V/3k ||l 4 (0,T;L 4 (ST)) 


2 C _i i 

< — 4 = 2 C/K 8 

Vk 


0 


as k —y 0 , where we used 3(7 + |) < 3 ^ for any 7 > 1 . 
Calculating P \, 

rT 


P = ~ 


f [ i’(t)pl(j) K {p K )ipl:X7v K dxdt 
Jo Jn 


/ 0 


V’(t)/ 0 ^(</>A'(p K ))V n : Vu K dxdt 


/ / l &{t)pZ<t>K(pK)<Pn : {V<I>k(Pk) ® u k ) dxdt 

Jo Jn 


= Pi 1 + P] 2 • 

We bound P 12 as follows 


1^121 = - 
r T 
1 0 


r 

0 -'O 


il>(t)pl</> K (pK)<Pn : < 8 > u k ) dxdt 

[ [ ^(t)pl{^(!>K{p K )) T p'n{v K )v K dxdt 

Jo Jn 


-7+3, 


<C7(n,M| L co)K s ||P/t ’ll L $ (0iT;L $ (n)) 

2(3* 1 1 

< -—^=K _ ^ = 2Ck5 -+ 0 
- \/iL 


N 4 ||V/?K || i 4( 0j 'r ; I, 4 (n)) ) II ( /4 '\/Pk||l 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


as ft ^ 0, where we used |</>* : (p K )- v /^l < ll^n( v fc) v fc||L°°(o.T ; i,°°(a)) < C(n) since 

is compactly supported, and 3(7 + 3 ) < 3 ^ for any 7 > 1 . 

Thanks to part b of Lemma 12.21 and 

¥>n(v K ) : Vu K = 4Vu K c^"(|v K | 2 )v K <g) v K + 2divu K ^(|v K | 2 ), 
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we have 


\Pn\ < 


+ 2 


and the term 


4 / [ ^(t)\(j) K \ 2 \p' y K \\Vu K \dxdt + 2 [ [ ip(t)\<fi K \ 2 \tp , n (\v K \ 2 )\\pJ.\\divu K \ dx dt 

Jo Jn Jo Jn 

f [ p K \Vu K \ 2 dxdt + C(\\ip\\ L °°) / f p 2-y_1 dx dt 
Jo Jn Jo Jn 

[ [ ^(tMK\ 2 \^n(W^\ 2 )\\pl\\^^\dxdt, 

Jo Jn 


< 4 I|V’IU°° 

,-T 


rT 


< 2 


^t)\<j)K\ 2 \v'n(\v K \ 2 )\\pl\\dwu K \dxdt 

[ [ K \ 2 )\p K \Il)u K \ 2 dx dt 

Jo Jn 


io Jn 

r -T 


+ 


C(M |l») [ f \M^k\ 2 )\pV 1 dx dt. 

Jo Jn 


Thus, 


|Ai|< 


1 2 dxdt 


The first right hand side term will be controlled by 


za° ^ (toluol 2 + + IV^I 2 - r 0 log_ po'j dx + 


2 Et 


(4.12) 


(4.13) 


/ / p K |Vu« 

Jo Jn 

+ 2 [ [ ^{t)\^ K \ 2 W n {W^\ 2 )\PnP^\ 2 dx dt (4.14) 

Jo Jn 

+ c (W\\l°°) [ [ (1 + <PniWK\ 2 ))p 2 K y ~ 1 dxdt: 

Jo Jn 


due to (11.161) : and the second right hand side term will be absorbed by the dispersion term 
Ai in (14.230 . By Lemma 14.11 we have 

[ [ (l + <Pn(l ^k\ 2 )pV^ dxdt^ f [ (l + ^' n (|u| 2 ))p 2l ~ l dxdt (4.15) 

Jo Jn Jo Jn 


as k —^ 0. 
Note that 


[ [ V’(*) < Pn( v «)( r o u « ; + rip K \u K \ 2 u K ) dx dt > 0, (4.16) 

Jo Jn 


so this term can be dropped directly. 
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We treat the other terms in F one by one, 

/ [ \^{t)(p' n (v K )p 2 K u K (l)'K(p K )divu K \ dxdt 

J o Jn 

i 

< C(n, '0)||PKU K || L 4((o i T;L 4 (n))llv / /’K^ VU K||L 2 ((O,T;L 2 (Q)) 

3 3 

x M'k(Pk)Vp^\\l°° Wpk llL 4 ((o,T ; L 4 (f 2 )) < C(n, V’)n® —> 0 
as n —> 0, where we used Sobolev inequality, and l^'xiP^yfPnl < H=; 

/ / |^(O^Uwc)PkV<)>a(/0 k )Buk| dxdt 

Jo Jn 

< C(n, Ip) '+MVP1 f K 4 \\VPk ||l 4 ((0,T;L4(Q))^ |Iv / Pk® )u k||l 2 (( 0 ,T;L 2 (Q))II/ 3 k ||l 4 (( 0 ,T;L 4 (Q)) 

K 4 V / 

< C(n,ip)K s -A 0 
as k o 0; 


(4.17) 


(4.18) 


c / / dxdt 

Jo Jn 

< 2 C(n,1p)n4 f« 4 ||Vp^ || L 4((0 i T;L 4 (r2)) N ) llv / NA v / p)T|| i 2 ( '( 0iT;L 2(Q)) ||/3« ||l 4 (0,T;Z, 4 (Q)) 


< 2C'(n,^)K 4 -a 0 

as k —^ 0, where we used |/9 k 0aO°k)I < 1- Finally 
K / |'0 / (i) ( / :, n(v« ; )(/)R(/5 K )V x /p):A v /^| dxdt 

Jo Jn 

< 2C'(n,V ; )K 4 3 11v||L4((o i r ; L 4 (n)) N ) WV^^-\ /r (F\\L 2 ((o,T-,L 2 (n)) ||p« ||l 4 (o,t ; l 4 (Q)) ^ 4 ' 20 ^ 


< 2 C(n,^)n^ -A 0 
as n —y 0. 


the term S K = 4>k(Pk)Pk.( Du k + = Si + S 2 , we calculate as follows 

T 

[ [ i/j{t)Si : V(^(v K )) dxdt 


1 0 Jn 

n 

1 0 Jn 

rT r 


^{t)4>K{p K )pJJ> u K : V(p' n (v K )) dx dt 


ip(t)[Vu K ip"(v K ,)p K ] : B>u K ((j) K (p K )) 2 dxdt 


Jo Jn 
= A 1 +^ 2 - 


■0 {t)p K <t>K(P k){ Uk (v k ))Du k V(0a- (p K ) ) dx dt 


(4.21) 
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For A \, by part a. of Lemma 12.21 we have 
rT 


= 


[ [ i/>(t)[Vu K ^(v K )/) K ] : Bu K ((j) K (p K )) 2 dxdt 

Jo Jn 

= 2 [ [ ^(t)if' n (\v K \ 2 )(^ K (p K )) 2 pJ])u K :Vu K dxdt 

Jo Jn 

+ 4/ / i/>(t) / 3 K (</>A'(/0 K )) 2 ^n(|v ft | 2 )(Vu K v K ® v K ) : Ou K dxdt 

Jo Jn 


(4.22) 


= An + A 12 . 


Notice that 


thus 


Bu k : Vu K = |Bu k | 2 , 


Ai > 2 


[ [ ip(t)<p' n (\v K \ 2 )((j) K (p K )) 2 p K \Bu K ,\ 2 dxdt 

Jo Jn 


/ / /Lt| Vu K 

Jo Jn 


(4.23) 


; dx dt, 


where we control A 12 
Ai 2 <4 [ T [ 


< 


0 J n 

Loo 


vr f f 

\il>{t)\r— r— ^p K \Vu K \ 2 dxdt < 4||^|U°° / / p K \ Vu K | 2 dx dt 

t + |v K | Jo Jn 

(X ( /5 °I U °I 2 + + I V \/Po| 2 “ r o lo g- , 0 ) dx + 2 E 0 ^j , 

thanks to (11.161) . For A 2 , thanks to (12.81) . we can control it as follows 

_ 3 

IA 2 1 < C(n, '0)||v / ^D u K||L2((o ) T;L 2 (a))l|PK||L 4 ((O,T;L 4 (Q)) 

b kVPkWl™ ((o,T)xa) 


X (k 4 II V/3k Hl4(( 0i T;L 4 (Q)))' 

c 


K 4 


(4.24) 


< 


Vk, 


— = Ck 8 -A 0 


K 4 


as ti —^ 0. Note that the first right hand side term of (14.231) has a positive sign and control 
the limit using from the pressure (14.141) . We need to treat the term related to S 2 , 

« [ [ V’(*)S2 : V(ip' n (v K ))dxdt 

Jo Jn 
rr ' _ 

Jc '^(*) Vu «^( v ' t ) : Vp^^k^Pk) 2 ky/Kdxdt ^ 25 ^ 

/ / V , (i)u K </ > x(pK)^n(v K )V p Ky /p^(f>' K (p K )Ay/p^dxdt 

Jo Jn 


= K 


+ K 
= .Bi + i?2 
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we control B± as follows 

l-Bil < Cin , iji ) ||-\/Pk||l 2 (o,t ; L 2 (fi)) II V^^\//^IIl 2 (o,t ; L 2 (n ))\/n 

< Cxh -a- o 


(4.26) 


as n —> 0, where we used W^kWl™ (o,T-.L°°{n)) < C. 
For B- 2 , we have 


\B 2 \ < C(n, , 0)k 4 ^4||Vp^|| L 4 (0iT;i 4 (n)) J 

X llP/t IIl 4 (0,T ;L 4 (fi)) II L 2 (0,T;L 2 (fi)) II (Pk)^k; Hi, 00 (0,T;L°°(S2)) 

< ChO -a- 0 


(4.27) 


as k —)■ 0. 

With (|4.6D - ()4.27I) . in particularly, letting n —> 0 in (13.2p . dropping the positive terms 
on the left side, we have the following inequality 


/ / ij/(t)pip n (u) dxdt 

Jo Jn 

(po|u 0 | 2 + + | V^Pol 2 - ro log_ Po^j dx 


+ 2 Eq 


+ 


^(0) [ pQ<p n (uo)dx + C(\\i/j\\ l °o) [ [ (1 + $ l (|u| 2 ))p 27 1 da;di, 

4o Jn 


which in turn gives us Lemma 14.21 


□ 


5. Recover the limits as n —>• oo. 

In this section, we aim at recovering the Mellet-Vasseur type inequality for the weak 
solutions at the approximation level of compressible Navier-Stokes equations by letting 
n —>• oo. In particular, we prove Theorem 11.11 by recovering the limit from Lemma 14.21 In 
this section, ( p , u) are the fixed weak solutions. 

Our task is to bound the right term of (14.51) . 

C'(IIV’IIl-) [ [ (l + £n(l u l 2 ))p 2l ~ X dxdt 
Jo Jn 

t / \ — / \- 

< C'(||^|| L oo) [ ( /(p 27_1_ 5)2^ dx ) ( [ p{l + ( p ' n i \ u \ 2 ))^ dx \ dt (5.1) 

Jo \Jn J \J n / 

rp 2^ « 

<C{ Wh-) [ ( [ ( P 2 ( [ p(2 + ln(l + |u| 2 )))ldx) 2 dt , 

Jo \Jn J \Jn J 
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where we used part c of Lemma 12.21 By (14.51) and (15.11) , we have 

- / ip'(t)pip n (u)dxdt < / p 0 p n (u 0 )dx 

Jo Jn Jn 

(/ ( p °l u °l 2 + + “ r ° log - /5 °) dx + 2 E o ) 


+ 8 U L» 


+ 


C(II^IU-) [ T ([ (p 2 ^- 1 --*) 

Jo \Jn 


2—5 

2 \ 2 
2—5 


p(2 + ln(l + |u | 2 )))« dx dt. 


Thanks to part e. of Lemma 12.21 and Monotone Convergence Theorem, we have 

— f f ’ip 1 (t)pip n (u) dx dt —> — f f ip'(t)p(l + |u| 2 ) ln(l + |u| 2 ) dx dt (5.2) 

Jo Jn Jo Jn 


as n —>• oo. 

Letting n —> oo, we have 

rT 


— f f ip'(t)p(l + | u | 2 ) ln(l + |ii| 2 ) dx dt < ip(0) f po(l + |u 0 | 2 )hi(l + |u 0 | 2 )da; 

Jo Jn Jn 

+ 8||i/’||l°° ^Po|u 0 | 2 + + iv^/p^l 2 - r 0 log_/9 0 ^ dx + 2£o^ (5.3) 


2-5 __ 

+ C f f [ (p 27-1- ^) 537 ^) ( [ p(2 + ln(l + |u| 2 ))« dx\ dt. 

Jo \Jn J \Jn J 


Taking 


ip(t) < 


= 1 

_ 1 _ uT 

2 e 
= 0 


if t < t — | 

if t— §<i<t+§ 

if t > t + f, 


then ()5.31) gives for every t > |, 
tea p(l + |u| 2 ) ln(l + |u| 2 ) dx^j dt 

< Po(l + |uo| 2 ) ln(l + |u 0 | 2 ) dx + 8 ^ ^Po|u 0 | 2 + + |V\/p^| 2 - r 0 log_ dx + 

5 

p(2 + ln(l + |u| 2 )) s dx J dt. 


2V 2 


+c io u, ( " 

This gives Theorem 11.11 thanks to the Lebesgue point Theorem. 

6. Recover the weak solutions 

The objective of this section is to apply Theorem 11.11 to prove Theorem 11.21 In partic¬ 
ular, we aim at establishing the existence of global weak solutions to (jl.ll) - (U.2j) by letting 
ro —>• 0 and r\ —» 0. Let r = tq = r\. we use (p r , u r ) to denote the weak solutions to (12.11) 
verifying Proposition 11.11 with k = 0. 
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By (11.71) and () 1.8 D . one obtains the following estimates, 

|| v / p^ u r || i oo( 0j T;L 2 (r2)) < C] 
|Pr||L°°(0,T;L 1 nL'i'(Q)) < Cj 
I \/P^V u r ||i 2 (o i T;L 2 (n)) < 

|v^:|| L oo (0,T;L 2 (f2)) < 


( 6 . 1 ) 


I^Pr^ 2 IIl 2 (0,T;L 2 (Q)) < C. 


Theorem ll.il gives us 


sup 


p r |u r p ln(l + |u r | 2 ) dx < C. (6.2) 

te[o,T] Jn 

It is necessary to remark that all above estimates on (16.11) and (16.21) are uniformly on r. 
Thus, we can make use of all estimates to recover the weak solutions by letting r —>• 0. 
Meanwhile, we have the following estimates from <tm 

f f r\u r \ 2 dx dt < C, 

(6.3) 

/ / rp r \u r \ A dx dt < C. 

Jo Jn 

To establish the existence of global weak solutions, we should pass to the limits as r —>• 0. 
Following the same line as in [32], we can show the convergence of the density and the 
pressure, prove the strong convergence of yfpju r in space Lf oc ((0, T) x fl), the convergence 
of the diffusion terms. We remark that Theorem o is the key tool to show the strong 
convergence of y/Jpu r . Here, we list all related convergence from [32]. In particular, 

yfpj —» yfp almost everywhere and strongly in Lf oc {((). T) x H)), 


Pr^p in C°(0,r;Lf oc (H)); 
the convergence of pressure 

pj —> p 7 strongly in L} oc ((0,T) X fi); 
the convergence of the momentum and y/Jpu r 

p r u r —>• pu strongly in L 2 (0,T; Lf oc (H)) for pE [1,3/2); 
yfpj\i r -A- y/pu strongly in L 2 oc ((0,T)xH); 
and the convergence of the diffusion terms 

/ 

p r Vu r —> pVu in 2) , 


(6.4) 


(6.5) 


( 6 . 6 ) 


p r \7 T u r -A- pV r u in 2) . 

It remains to prove that terms ru r and r/? r |u r | 2 u r tend to zero as r 
test function, then we estimate the term ru r 


(6.7) 

0. Let il* be any 


ff 


rupt/j dx dt 


< 


f T f I!, 

/ / r 2 r 2 1 u, r 

7o Jn 


dx dt 


< \A’||\A : U r || Z/ 2((o i T)xn)IIV ; llL 2 ((0,T)xQ) 0 


( 6 . 8 ) 
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as r —>• 0, due to (16.31) . 

We also estimate rp r |u r | 2 u r as follows 

/ / rp r | u r | 2 u r t/; dx dt 

Jo Jn 


(6.9) 


< \/r\\y/r u r | ^ IIl 2 ((0,T)xO) II VPr W - ||l°° (0,T ;L 2 (f2)) IIV 1 III. 00 ((0,T)x fi) “> 0 

as r —> 0. 

The global weak solutions to (12.11) verifying ProDosition ll.il with n = 0 is in the following 
sense, that is, (p r ,u r ) satisfy the following weak formulation 


/ p r \i r ■ %p dx\\J^ — / / p r \i r ipt dx dt — / / p r u r (8) u r : Vi/’ dx dt 

Jn Jo Jn Jo Jn 

— / pjdiv'if; dx dt — 

Jo Jn 

rT 


n 

'o Jn 


= —r 


fu 


ip dxdt — r 


pBu r : V tp dxdt 
/ / p r \u. r \ 2 u r ^ dx dt, 

Jo Jn 


( 6 . 10 ) 


where ^p is any test function. 

Letting r —>• 0 in the weak formulation (16.101) . and applying (I6.4D - (I6.9I) . one obtains that 


pu -ipdx |t = o - 


n 

Jo Jn 


puipt dx dt — 


to Jn 


p^divip dx dt — 


o Jn 

f f pDu : Vtp dx dt = 0. 
Jo Jn 


pu <g) u : V'L dx dt 


( 6 . 11 ) 


Thus we proved Theorem ll.2l for any initial value (poi u o) verifying (11.51) with the additional 
condition po > This last condition can be dropped using |[32| . 
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